Effects of Donor Ages and Propagation Methods on Seedling Growth of Platycladus orientalis (L.) Franco in Winter.

To evaluate the effects of donor ages on growth and stress resistance of 6-year-old seedlings propagated from 5-, 2000-, and 3000-year-old Platycladus orientalis donors with grafting, cutting, and seed sowing, growth indicators and physiological and transcriptomic analyses were performed in 6-year-old seedlings in winter. Results showed that basal stem diameters and plant heights of seedlings of the three propagation methods decreased with the age of the donors, and the sown seedlings were the thickest and tallest. The contents of soluble sugar, chlorophyll, and free fatty acid in apical leaves of the three propagation methods were negatively correlated with donor ages in winter, while the opposite was true for flavonoid and total phenolic. The contents of flavonoid, total phenolic, and free fatty acid in cutting seedlings were highest in the seedlings propagated in the three methods in winter. KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of differentially expressed genes showed phenylpropanoid biosynthesis and fatty acid metabolism pathways, and their expression levels were up-regulated in apical leaves from 6-year-old seedlings propagated from 3000-year-old P. orientalis donors. In addition, hub genes analysis presented that C4H, OMT1, CCR2, PAL, PRX52, ACP1, AtPDAT2, and FAD3 were up-regulated in cutting seedlings, and the gene expression levels decreased in seedlings propagated from 2000- and 3000-year-old donors. These findings demonstrate the resistance stability of cuttings of P. orientalis and provide insights into the regulatory mechanisms of seedlings of P. orientalis propagated from donors at different ages in different propagation methods against low-temperature stress.

Seeing through the hedge: Phylogenomics of Thuja (Cupressaceae) reveals prominent incomplete lineage sorting and ancient introgression for Tertiary relict flora.

The Eastern Asia (EA) - North America (NA) disjunction is a well-known biogeographic pattern of the Tertiary relict flora; however, few studies have investigated the evolutionary history of this disjunction using a phylogenomic approach. Here, we used 2369 single copy nuclear genes and nearly full plastomes to reconstruct the evolutionary history of the small Tertiary relict genus Thuja, which consists of five disjunctly distributed species. The nuclear species tree strongly supported an EA clade Thuja standishii-Thuja sutchuenensis and a "disjunct clade", where western NA species T. plicata is sister to an EA-eastern NA disjunct Thuja occidentalis-Thuja koraiensis group. Our results suggested that the observed topological discordance among the gene trees as well as the cytonuclear discordance is mainly due to incomplete lineage sorting, probably facilitated by the fast diversification of Thuja around the Early Miocene and the large effective population sizes of ancestral lineages. Furthermore, approximately 20% of the T. sutchuenensis nuclear genome is derived from an unknown ancestral lineage of Thuja, which might explain the close resemblance of its cone morphology to that of an ancient fossil species. Overall, our study demonstrates that single genes may not resolve interspecific relationships for disjunct taxa, and that more reliable results will come from hundreds or thousands of loci, revealing a more complex evolutionary history. This will steadily improve our understanding of their origin and evolution.

Plastome Phylogenomic and Biogeographical Study on Thuja (Cupressaceae).

Investigating the biogeographical disjunction of East Asian and North American flora is key to understanding the formation and dynamics of biodiversity in the Northern Hemisphere. The small Cupressaceae genus Thuja, comprising five species, exhibits a typical disjunct distribution in East Asia and North America. Owing to obscure relationships, the biogeographical history of the genus remains controversial. Here, complete plastomes were employed to investigate the plastome evolution, phylogenetic relationships, and biogeographic history of Thuja. All plastomes of Thuja share the same gene content arranged in the same order. The loss of an IR was evident in all Thuja plastomes, and the B-arrangement as previously recognized was detected. Phylogenomic analyses resolved two sister pairs, T. standishii-T. koraiensis and T. occidentalis-T. sutchuenensis, with T. plicata sister to T. occidentalis-T. sutchuenensis. Molecular dating and biogeographic results suggest the diversification of Thuja occurred in the Middle Miocene, and the ancestral area of extant species was located in northern East Asia. Incorporating the fossil record, we inferred that Thuja likely originated from the high-latitude areas of North America in the Paleocene with a second diversification center in northern East Asia. The current geographical distribution of Thuja was likely shaped by dispersal events attributed to the Bering Land Bridge in the Miocene and subsequent vicariance events accompanying climate cooling. The potential effect of extinction may have profound influence on the biogeographical history of Thuja.

Development of microsatellite markers for the Japanese endemic conifer Thuja standishii and transfer to other East Asian species.

OBJECTIVE: Design polymorphic microsatellite loci that will be useful for studies of the genetic diversity, gene-flow and reproduction in the Japanese endemic conifer Thuja standishii and test the transferability of these loci to the two other East Asian species, T. sutchuenensis and T. koraiensis. RESULTS: Fifteen loci were developed which displayed 3 to 21 alleles per locus (average = 9.2) among 97 samples from three populations of T. standishii. Observed heterozygosity for all samples varied between 0.33 and 0.75 (average = 0.54) while expected heterozygosity values were higher with an average over the 15 loci of 0.62 (0.37-0.91). Low multi-locus probability of identity values (< 0.00002) indicate that these markers will be effective for identifying individuals derived from clonal reproduction. All 15 loci amplified in 13 samples of T. sutchuenensis, the sister species of T. standishii, with 1 to 11 alleles per locus (average = 4.33) while 13 loci amplified in four samples of the more distantly related T. koraiensis with 1 to 5 alleles per locus (average = 2.15).

Exploring the Phylogeography of Ancient Platycladus orientalis in China by Specific-Locus Amplified Fragment Sequencing.

Platycladus orientalis (i.e., Chinese thuja) is famous for its lifespan spanning hundreds, and even thousands, of years. Most ancient P. orientalis populations are widely distributed in China, with accessible historical records, making them valuable genetic resources. In this study, the distribution pattern of ancient P. orientalis in China was analyzed based on 13 bioclimatic factors. Additionally, a specific-locus amplified fragment (SLAF) sequencing method was applied to detect single nucleotide polymorphisms (SNPs) in the genomes of 100 accessions from 13 populations. The resulting data revealed that the suitable areas for the distribution of ancient P. orientalis populations were accurately predicted with four main climatic factors. A total of 81,722 SNPs were identified from 461,867 SLAFs for 100 individuals, with an average sequencing depth of 10.11-fold and a Q30 value of 82.75%. The pair-wise genetic distance and genetic differentiation of 13 populations indicated that the BT-T population exhibited the largest divergence from the other populations. A neighbor-joining phylogenetic tree suggested the relationship between many individuals was inconsistent with the geographical location, possibly indicative of a history of transplantation and cultivation. All individuals were clustered into nine genotypes according to a structural analysis and the relationships between individuals were clarified in phylogenetic trees. This study highlights the importance of the de novo genome sequencing of ancient P. orientalis and may provide the basis for the conservation of P. orientalis genetic resources, the identification of supergene families, and the evaluation of related genetic resources.

The gymnosperm cytochrome P450 CYP750B1 catalyzes stereospecific monoterpene hydroxylation of (+)-sabinene in thujone biosynthesis in western redcedar.

Western redcedar (WRC; Thuja plicata) produces high amounts of oxygenated thujone monoterpenoids associated with resistance against herbivore feeding, particularly ungulate browsing. Thujones and other monoterpenoids accumulate in glandular structures in the foliage of WRC. Thujones are produced from (+)-sabinene by sabinol and sabinone. Using metabolite analysis, enzyme assays with WRC tissue extracts, cloning, and functional characterization of cytochrome P450 monooxygenases, we established that trans-sabin-3-ol but not cis-sabin-3-ol is the intermediate in thujone biosynthesis in WRC. Based on transcriptome analysis, full-length complementary DNA cloning, and characterization of expressed P450 proteins, we identified CYP750B1 and CYP76AA25 as the enzymes that catalyze the hydroxylation of (+)-sabinene to trans-sabin-3-ol. Gene-specific transcript analysis in contrasting WRC genotypes producing high and low amounts of monoterpenoids, including a glandless low-terpenoid clone, as well as assays for substrate specificity supported a biological role of CYP750B1 in α- and ß-thujone biosynthesis. This P450 belongs to the apparently gymnosperm-specific CYP750 family and is, to our knowledge, the first member of this family to be functionally characterized. In contrast, CYP76AA25 has a broader substrate spectrum, also converting the sesquiterpene farnesene and the herbicide isoproturon, and its transcript profiles are not well correlated with thujone accumulation.

Genetic diversity of the critically endangered Thuja sutchuenensis revealed by ISSR markers and the implications for conservation.

Thuja sutchuenensis Franch. is a critically endangered plant endemic to the North-East Chongqing, China. Genetic variation was studied to assess the distribution of genetic diversity within and among seven populations from the single remnant locations, using inter-simple sequence repeat (ISSR) markers. A total of 15 primers generated 310 well defined bands, with an average of 20.7 bands per primer. The seven populations revealed a relatively high level of genetic diversity in the species. The percentage of polymorphic bands, Nei's gene diversity and Shannon's information index at the population and species level were 76.1%, 0.155, 0.252 and 100%, 0.165, 0.295, respectively. A low level of genetic differentiation among populations (G(ST) = 0.102), in line with the results of Analyses of Molecular Variance (AMOVA), and a high level of gene flow (N(m) = 4.407) were observed. Both the Unweighted Pair Group Method with Arithmatic Mean (UPGMA) cluster analysis and Principal Coordinates Analysis (PCoA) supported the grouping of all seven populations into two groups. In addition, Mantel test revealed no significant correlation between genetic and geographical distances (r = 0.329, p = 0.100). The low genetic differentiation among populations implies that the conservation efforts should aim to preserve all the extant populations of this endangered species.

Identification of genes in Thuja plicata foliar terpenoid defenses.

Thuja plicata (western redcedar) is a long-lived conifer species whose foliage is rarely affected by disease or insect pests, but can be severely damaged by ungulate browsing. Deterrence to browsing correlates with high foliar levels of terpenoids, in particular the monoterpenoid α-thujone. Here, we set out to identify genes whose products may be involved in the production of α-thujone and other terpenoids in this species. First, we generated a foliar transcriptome database from which to draw candidate genes. Second, we mapped the storage of thujones and other terpenoids to foliar glands. Third, we used global expression profiling to identify more than 600 genes that are expressed at high levels in foliage with glands, but can either not be detected or are expressed at low levels in a natural variant lacking foliar glands. Fourth, we used in situ RNA hybridization to map the expression of a putative monoterpene synthase to the epithelium of glands and used enzyme assays with recombinant protein of the same gene to show that it produces sabinene, the monoterpene precursor of α-thujone. Finally, we identified candidate genes with predicted enzymatic functions for the conversion of sabinene to α-thujone. Taken together, this approach generated both general resources and detailed functional characterization in the identification of genes of foliar terpenoid biosynthesis in T. plicata.

Genetic diversity and differentiation of core vs. peripheral populations of eastern white cedar, Thuja occidentalis (Cupressaceae).

PREMISE OF THE STUDY: Geographically peripheral (marginal) populations are expected to have lower genetic diversity and higher genetic differentiation than geographically core (central) populations as a result of supposedly lower effective population size (N(e)) and higher genetic drift, founder effect, fragmentation, and isolation in peripheral than in core populations. Here we address this issue for a long-lived plant species, eastern white cedar (Thuja occidentalis). METHODS: Genetic diversity and population structure of 13 natural populations of eastern white cedar from its Canadian eastern peripheral and core natural ranges in New Brunswick, Nova Scotia, and Prince Edward Island were studied using six nuclear microsatellite DNA markers. KEY RESULTS: The core populations of eastern white cedar had significantly higher allelic diversity (mean A = 8.83, A(r) = 8.13, A(e) = 4.03) and N(e) (428) than the peripheral populations (A = 6.64, A(r) = 6.15, A(e) = 3.12, N(e) = 198). However, expected heterozygosity was similar in the core (H(e) = 0.64) and peripheral (H(e) = 0.60) populations. Genetic differentiation was significantly higher among the peripheral (F(ST) = 0.089) than among the core (F(ST) = 0.032) populations. No genetic differentiation (F(ST)/Φ(RT) = 0.000) was detected between core and peripheral regions. CONCLUSIONS: Peripheral populations have significantly lower N(e) and genetic diversity in terms of allelic diversity (richness) and significantly higher genetic differentiation than the core populations of eastern white cedar in its Canadian eastern range. However, core and peripheral populations have similar levels of expected heterozygosity. Implications for conservation of eastern white cedar genetic resources are discussed.

Higher fine-scale genetic structure in peripheral than in core populations of a long-lived and mixed-mating conifer--eastern white cedar (Thuja occidentalis L.).

BACKGROUND: Fine-scale or spatial genetic structure (SGS) is one of the key genetic characteristics of plant populations. Several evolutionary and ecological processes and population characteristics influence the level of SGS within plant populations. Higher fine-scale genetic structure may be expected in peripheral than core populations of long-lived forest trees, owing to the differences in the magnitude of operating evolutionary and ecological forces such as gene flow, genetic drift, effective population size and founder effects. We addressed this question using eastern white cedar (Thuja occidentalis) as a model species for declining to endangered long-lived tree species with mixed-mating system. RESULTS: We determined the SGS in two core and two peripheral populations of eastern white cedar from its Maritime Canadian eastern range using six nuclear microsatellite DNA markers. Significant SGS ranging from 15 m to 75 m distance classes was observed in the four studied populations. An analysis of combined four populations revealed significant positive SGS up to the 45 m distance class. The mean positive significant SGS observed in the peripheral populations was up to six times (up to 90 m) of that observed in the core populations (15 m). Spatial autocorrelation coefficients and correlograms of single and sub-sets of populations were statistically significant. The extent of within-population SGS was significantly negatively correlated with all genetic diversity parameters. Significant heterogeneity of within-population SGS was observed for 0-15 m and 61-90 m between core and peripheral populations. Average Sp, and gene flow distances were higher in peripheral (Sp = 0.023, σg = 135 m) than in core (Sp = 0.014, σg = 109 m) populations. However, the mean neighborhood size was higher in the core (Nb = 82) than in the peripheral (Nb = 48) populations. CONCLUSION: Eastern white cedar populations have significant fine-scale genetic structure at short distances. Peripheral populations have several-folds higher within-population fine-scale genetic structure than core populations. Anthropogenic disturbances and population fragmentation presumably have significant effects on fine-scale genetic structure in eastern white cedar. Core populations have higher neighborhood size than peripheral populations, whereas gene flow distances are higher in peripheral than in core populations. The results of our study contribute to the knowledge of poorly-understood spatial genetic structure of core versus peripheral populations in plants. As well, the information is of significance for conservation of genetic resources of eastern white cedar and perhaps of other long-lived forest trees with mixed-mating system.

Activity, polypeptide and gene identification of thylakoid Ndh complex in trees: potential physiological relevance of fluorescence assays.

Three evergreen (Laurus nobilis, Viburnum tinus and Thuja plicata) and two autumnal abscission deciduous trees (Cydonia oblonga and Prunus domestica) have been investigated for the presence (zymogram and immunodetection) and functionality (post-illumination chlorophyll fluorescence) of the thylakoid Ndh complex. The presence of encoding ndh genes has also been investigated in T. plicata. Western assays allowed tentative identification of zymogram NADH dehydrogenase bands corresponding to the Ndh complex after native electrophoresis of solubilized fractions from L. nobilis, V. tinus, C. oblonga and P. domestica leaves, but not in those of T. plicata. However, Ndh subunits were detected after SDS-PAGE of thylakoid solubilized proteins of T. plicata. The leaves of the five plants showed the post-illumination chlorophyll fluorescence increase dependent on the presence of active Ndh complex. The fluorescence increase was higher in autumn in deciduous, but not in evergreen trees, which suggests that the thylakoid Ndh complex could be involved in autumnal leaf senescence. Two ndhB genes were sequenced from T. plicata that differ at the 350 bp 3' end sequence. Comparison with the mRNA revealed that ndhB genes have a 707-bp type II intron between exons 1 (723 bp) and 2 (729 bp) and that the UCA 259th codon is edited to UUA in mRNA. Phylogenetically, the ndhB genes of T. plicata group close to those of Metasequoia, Cryptomeria, Taxodium, Juniperus and Widdringtonia in the cupresaceae branch and are 5' end shortened by 18 codons with respect to that of angiosperms.

Relationship of dirigent protein and 18s RNA transcript localization to heartwood formation in western red cedar.

Western red cedar (Thuja plicata) heartwood contains abundant amounts of structurally complex plicatic acid-derived lignans that help confer protective properties and longevity to this tissue type. Although the lignan biochemical entry point is dirigent protein-mediated, the formation of heartwood and its associated lignans in some species remains poorly understood due to technical difficulties of working with the former. To begin to address such questions, this study therefore focused on the anatomical localization of dirigent protein and 18s rRNA (control) gene transcripts within recalcitrant woody tissues, including heartwood. This in situ mRNA hybridization approach enabled detection of dirigent protein transcripts in cork cambia, vascular cambia and ray parenchyma cells of the sapwood, but not the heartwood under the conditions employed. By contrast, the hybridization of the 18s rRNA (control) transcript resulted in its detection in all tissue types, including radial parenchyma cells of apparently preformed heartwood. Application of in situ hybridization to such recalcitrant tissues thus demonstrates the utility of this technique in identifying specific cell types involved in heartwood formation, as well as the relationship of dirigent protein localization to that of heartwood metabolite generation.

Reticulate evolution in Thuja inferred from multiple gene sequences: implications for the study of biogeographical disjunction between eastern Asia and North America.

The eastern Asia-North America disjunction is one of the most interesting biogeographical patterns, but its formation is still in much debate. Here nucleotide sequences of five cpDNA regions, nrDNA ITS and two low-copy nuclear genes (LEAFY, 4CL) were employed to reconstruct the phylogeny and to explore the historical biogeography of Thuja, a typical eastern Asia-North America disjunct genus. High topological discordance was observed between chloroplast and nuclear gene trees, even between different nuclear gene trees, suggesting that Thuja could have a reticulate evolutionary history due to multiple interspecific hybridization events. The eastern Asian species Thuja koraiensis might have obtained its chloroplast genome from the eastern North American species T. occidentalis by chloroplast capture, while the western North American species T. plicata is very likely to have inherited a recombinant cpDNA. Based on the phylogenetic analysis of multiple genes, DIVA-reconstruction of the distribution history, molecular clock estimation and fossil data, we inferred that Thuja could have originated from the high-latitude areas of North America in the Paleocene or earlier with subsequent expansion into eastern Asia through the Bering Land Bridge. The two eastern Asia species T. standishii and T. sutchuenensis have a sister relationship, and their split could have occurred in the Oligocene or early Miocene. In the present study, the selection of molecular markers in biogeographic studies was also discussed. Since most previous studies on the eastern Asia and North America disjunction are based on uniparentally inherited cpDNA and (or) directly sequenced nrDNA ITS data, the historical reticulate evolution in the studied groups might have been underestimated. Therefore, we suggest that multiple genes from different genomes, especially low-copy nuclear genes, be used in this research area in the future.

Somatic mutations at microsatellite loci in western Redcedar (Thuja plicata: Cupressaceae).

A per-generation somatic mutation rate for microsatellites was estimated in western redcedar (Thuja plicata, Donn ex D. Don.: Cupressaceae). A total of 80 trees representative of the average size and age of reproductive trees were sampled in four natural populations in southwestern British Columbia. Samples of bulked haploid megagametophytes were collected from two or three positions on each tree, assuming that the collections were far enough apart that the same mutant sector was not sampled twice. All samples were genotyped at eight microsatellite loci. A single mutation corresponding to a stepwise increase in one dinucleotide repeat was detected. The estimated mutation rate for microsatellites was 6.3 x 10(-4) mutations per locus per generation (or 3.1 x 10(-4) per allele per generation), with a 95% confidence interval of 3.0 x 10(-5) to 4.0 x 10(-3) mutations per locus. Somatic mutations can contribute to a greater mutational load in trees, as compared to shorter lived plants, and genotypic mosaics within an individual have important implications for plant defense strategies and plant evolution.

Monolignol radical-radical coupling networks in western red cedar and Arabidopsis and their evolutionary implications.

The discovery of a nine-member multigene dirigent family involved in control of monolignol radical-radical coupling in the ancient gymnosperm, western red cedar, suggested that a complex multidimensional network had evolved to regulate such processes in vascular plants. Accordingly, in this study, the corresponding promoter regions for each dirigent multigene member were obtained by genome-walking, with Arabidopsis being subsequently transformed to express each promoter fused to the beta-glucuronidase (GUS) reporter gene. It was found that each component gene of the proposed network is apparently differentially expressed in individual tissues, organs and cells at all stages of plant growth and development. The data so obtained thus further support the hypothesis that a sophisticated monolignol radical-radical coupling network exists in plants which has been highly conserved throughout vascular plant evolution.

The western red cedar (Thuja plicata) 8-8' DIRIGENT family displays diverse expression patterns and conserved monolignol coupling specificity.

The isolation and characterization of a multigene family of the first class of dirigent proteins (namely that mainly involved in 8-8' coupling leading to (+)-pinoresinol in this case) is reported, this comprising of nine western red cedar (Thuja plicata) DIRIGENT genes (DIR1-9) of 72-99.5% identity to each other. Their corresponding cDNA clones had coding regions for 180-183 amino acids with each having a predicted molecular mass of ca. 20 kDa including the signal peptide. Real time-PCR established that the DIRIGENT isovariants were differentially expressed during growth and development of T. plicata (P < 0.05). The phylogenetic relationships and the rates and patterns of nucleotide substitution suggest that the DIRIGENT gene may have evolved via paralogous expansion at an early stage of vascular plant diversification. Thereafter, western red cedar paralogues have maintained an high homogeneity presumably via a concerted evolutionary mode. This, in turn, is assumed to be the driving force for the differential formation of 8-8'-linked pinoresinol derived (poly)lignans in the needles, stems, bark and branches, as well as for massive accumulation of 8-8'-linked plicatic acid-derived (poly)lignans in heartwood.

Linking deer browsing and terpene production among genetic identities in Chamaecyparis nootkatensis and Thuja plicata (Cupressaceae).

To investigate whether differential herbivore browsing reflects genetic variation in plant defense expression, variation in needle terpenes and damage caused by black-tailed deer (Odocoileus hemionus) was analyzed on yellow-cedar (Chamaecyparis nootkatensis) and western redcedar (Thuja plicata). In a 100-genet yellow-cedar population, three genets that were heavily browsed and had extremely low levels of monoterpenes (0-0.36% dry matter), sesquiterpenes, and diterpenes were compared to unbrowsed genets (0.85-3.83% monoterpenes in dry matter). These differences were maintained in individuals protected from browsing, suggesting genetically based variation in constitutive terpene production. In western redcedar, heavily browsed trees had significantly lower total monoterpene concentrations (1.69% dry matter) than lightly browsed trees (3.32% dry matter). One heavily browsed tree expressed no monoterpenes. No differences were found for diterpenes. In both species, the genotypes with extremely low monoterpene concentrations came from the same open-pollinated families.